Adiabatic nonlinear probes of one-dimensional Bose gases 
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We discuss two complimentary problems: adiabatic loading of one-dimensional bosons into an 
optical lattice and merging two one-dimensional Bose systems. Both problems can be mapped 
to the sine-Gordon model. This mapping allows us to find power-law scalings for the number of 
excitations with the ramping rate in the regime where the conventional linear response approach 
fails. We show that the exponent of this power law is sensitive to the interaction strength. In 
particular, the response is larger, or less adiabatic, for strongly (weakly) interacting bosons for the 
loading (merging) problem. Our results illustrate that in general the nonlinear response to slow 
relevant perturbations can be a powerful tool for characterizing properties of interacting systems. 



Dynamics of one-dimensional (ID) cold gases 0,0,0] is 
characterized by a strong response to external perturba- 
tions [1, 01 compared to higher-dimensional systems, due 
to the large density of low-energy states. For example, 
for slow processes the excitations are strongly enhanced 
and can qualitatively affect the dynamics and lead to the 
breakdown of adiabatic approximation [||. Understand- 
ing slow dynamics has two-fold implications. First, such 
dynamics can be used as a probe of interacting systems in 
the regimes where the linear response (Kubo-type) anal- 
ysis fails. Second, this understanding can help to devise 
optimal ways for minimizing heating in the system during 
slow processes. This minimization can be important for 
various applications including adiabatic preparation of 
strongly correlated systems and adiabatic quantum com- 
puting. In this Letter we will demonstrate both such 
implications using two related problems of loading ID 
interacting bosons into an optical lattice and merging 
(or splitting) two ID condensates. 

The physical parameter that defines the properties of 
ID bosonic systems is the ratio of the interaction and 
kinetic energies 7 = mc/(h 2 p), where m is the mass 
of a particle, c is the interaction strength Q, and p is 
the density of particles. In the spatially uniform sys- 
tem, the energy spectrum is gaplcss in the whole range 
of the parameter 7, which at low energies is reflected in 
the Luttinger liquid description [8[ . In the extreme limit 
of strong interactions, 7 ^> 1, Tonks- Girardeau (TG) 
gas 0, the particles behave as hard-core spheres, with 
the repulsion playing a role similar to the Pauli exclusion 
principle. This regime has been recently realized in cold 
atoms 0, El. 

External perturbations can significantly modify the 
phase structure of bosonic gases. A bosonic sys- 
tem undergoes the transition from superfluid to Mott- 
insulator 12, [l3j], driven by the competition of the ki- 
netic and interaction energy: in deep lattices the insu- 
lating phase is realized when the interaction energy be- 
comes larger than the tunneling (kinetic) energy. This 
mean-field picture of the transition is modified in one 
dimension. If the repulsion between particles is strong 
enough, an arbitrarily small periodic potential com- 



FIG. 1: Two processes equivalently described by the sine- 
Gordon model: a) loading of a ID Bose gas into a commen- 
surate optical lattice, b) merging of two ID Bose gases. 



mensurate with the boson density stabilizes the Mott- 
insulator phase 



14|, |15|. At weaker interactions, the per- 



turbation is irrelevant, and the system is gapless. 

We consider the response of the system to a linear ramp 
of the lattice amplitude V(t) with time, as it is sketched 
in Fig. [IJi) . A similar problem has been addressed experi- 
mentally in the case of a three-dimensional Bose-Einstein 
condensate [HI . The main result of our work is the scal- 
ing relation between the total number of excitations and 
the ramping rate S = \dtV\: 



5*1, 77= 1/(3- K), 0<K<2, 



(1) 



where K is the Luttinger parameter. For repulsive bosons 
K ?s 1 + 4/7 in the strongly interacting case, 7 3> 1, and 
K w 7r/^/7 for weak interactions, 7 < 1 [13, EH- Tnc 
scaling ([I]) is valid provided that the loading time is much 
longer than the transition time at which the excitations 
become effectively frozen (see Eq.© below). The scal- 
ing exponent r\ increases towards K = 2, where it sat- 
urates to unity. The critical point K = 2 separates the 
gapped phase at K < 2 where a weak lattice is relevant 
from the gapless phase at K > 2 where it is irrelevant. 
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FIG. 2: Scaling exponent r\ for the excitations n cx ~ S v as 
a function of the Luttinger parameter K. The Kosterlitz- 
Thouless (K-T) critical point at K = 2 separates gapped and 
gapless phases. The dashed line at K > 2 indicates linear 
response to loading of the lattice. The behavior of r\ in the 
vicinity of K — 2 is a sketch: a rounding-off of the cusp 
is expected for a finite system. Inset: The total number of 
excitations as a function of the ramping rate S for K = 1 
(dashed red line) and K < 2 (solid blue line). 



The response of the system to the perturbation is not 
universal in the latter case since it depends on high en- 
ergy cut-off of the effective Luttinger liquid description. 
Nonetheless, for a finite system, from the linear response 
approach one expects rj = 1 (indicated by the dashed line 
in Fig. \2§ ■ The relation (JTJ) can be generalized to an arbi- 
trary power r of the ramp of the lattice potential in time: 
n cx oc fi r /( 2+r - K \ For the linear ramp r = 1 this result 
reduces to Eq. ([1]) and for the exponential ramp r — > oo 
we obtain linear scaling as generally expected 1|| 2^| in 
ID systems with the linear spectrum. 

The scaling relation |T]) is obtained by mapping our 
problem to the sine-Gordon (SG) model [2l[ with time- 
dependent coupling. The same formalism describes a 
different system of two identical coupled ID Bose sys- 
tems (see Fig. [TJj)). The tunneling of particles leads 
to the Josephson-like energy term depending on the rela- 
tive phase of the condensates [22| , and this system can be 
also mapped to the sine-Gordon model. Merging the con- 
densates by adiabatic increase of the coupling produces 
defects in the relative phase. This process is character- 
ized by the scaling relation similar to |T]) with the sub- 
stitution K — * 1/(2K). Therefore, contrary to the prob- 
lem of loading into a lattice, the response of the system 
is more adiabatic for strongly-interacting bosons. This 
result is intuitively clear since stronger interactions be- 
tween bosons lead to larger phase fluctuations and there- 
fore to weaker effect of the tunneling 23|, l24[ . Another 
system where our results (TT|) apply is a ID Bose gas with 



long-range interactions. There the transition between the 
Mott-insulator and the Haldane-insulator phases is also 
described by the sine-Gordon model [IH with the tuning 
parameter being equivalent to the strength of the peri- 
odic potential in the loading problem. 

Let us focus on a specific problem of loading ID bosons 
into an optical lattice. At low energies the system with- 
out the lattice can be described within the Luttinger liq- 
uid formalism [§, The optical lattice, Vt(a;,t) = 
Vo(t) cos(27rx/a), commensurate with the boson density, 
causes interference between particles moving in opposite 
directions. These processes generate the well-known sine- 
Gordon Hamiltonian 

1 

2 

where the length and the time are measured in units of 
the lattice spacing a and a/v s , respectively, and the en- 
ergies are rescaled by a factor hv s /a, with the sound 
velocity v s = vf/K proportional to the Fermi velocity 
vp = tmp/m. The conjugate variables $(a;) and 6{x) 
describe the fluctuations of the supcrfluid phase and the 
density, respectively fhsj . The only important parameter 
that characterizes the sine-Gordon model is = 2\J-kK: 
the gapped phase corresponds to K < 2 and the gapless 



21|: 



n = - dx [<s>{xf + {d x ef - w(t) cos(pe)] , (2) 



does to K > 2 [2l|. The spectrum of the gapped phase 



consists of solitons and anti-solitons in the so-called re- 
pulsive regime at 1 < K < 2, and also their bound states, 
breathers, in the attractive regime at < K < 1 (2||. 
Solitons (anti-solitons) have topological charge +1 (— 1), 
while breathers have zero topological charge. For repul- 
sive bosons (K > 1) it is possible to cross the quantum 
critical point at K = 2 upon changing the interaction 
strength between particles. We note that in a gapless 
regime, there is a critical depth of the lattice at which 
a transition to a gapped Mott-insulator state occurs de- 
scribed by the effective renormalized K = 2 [3]. Other 
models, which we mentioned earlier correspond to other 
domains of K or equivalently (3. 

We assume a linear time-dependence of the lattice 
amplitude V{t) = VfOit, where Vf is the final value 
of the lattice, so that the loading time is restricted to 
< t < tf = 1/a. In the gapped regime the excitations 
are mainly created at initial stage, with their number be- 
ing effectively unchanged after some characteristic time 
scale related to the ramping rate S = \Vfa\ (see below). If 
the loading time is larger than this scale then our results 
do not depend on tf and we are able to make universal 
predictions. 

To study the excitation probability we employ the per- 
turbation theory in the ramping rate [27| : 
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E 



dt(p,t\d t \Q,t)e i ti dt '\. u * ( - t ' ) - Uo{ - t ' 



(3) 



where we assumed that the system is initially in its 
ground state, and the upper limit of integration over the 
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loading time was set to infinity. Here u> p (t) — uio(t) is 
the excitation energy defined in the instantaneous basis 
\p,t) with respect to the ground state \0,t). The fact 
that the matrix element of the perturbation appears in 
the instantaneous basis highlights the difference of this 
approach with the conventional linear response, in which 
unperturbed basis is employed. One can show that in the 
regimes where our approach gives nontrivial scaling, the 
linear response in the coupling V diverges. 

The time-dependent perturbation in Eq. (|3|) conserves 
the topological charge and the total momentum. There- 
fore in the adiabatic process the only relevant excitations 
are the soliton-antisoliton pairs created from the vacuum 
state. The instantaneous eigen-state \p, t) in Eq. ([4]) cor- 
responds to a pair with momenta (p, —p) , and the energy 
2 with a gap equal to twice the soliton 



2vV 



mi 



p 

mass m s . Therefore the matrix element (p,t\dt\0,t) for 
the time-dependent perturbation can be explicitly writ- 
ten as follows: 



(p,t\dt\0,t)=AdtV 



(0, t\A± {p)A T (-p) cos P6\0, t) 
u (t) - w p (t) 



(4) 



where A\_(p) is the operator that creates a soliton (an- 
tisoliton) state with momentum p. The matrix element 
appearing in Eq. (fj| is called the form- factor of cos(/30), 
and it is known from the literature (see Ref. [28^. 

The scaling behavior of n ex (^) is found from the anal- 
ysis of the phase factor in Eq.([3]). The dependence of the 
soliton mass m s ~ yVC 2-1 *") on the loading amplitude 
V extracted from the exact result of Ref. [29|, suggests 
rescaling the variables: 



(5) 



which leads to the invariant form of the phase 
/ dt' [u> p (t') — uio(t')] = ip(q,r). Similarly one obtains 
dt(p,t\d t \0,t) = giqr^^-^dV/V, where the function 
g can be expressed through the form- factors of Ref. [28| . 

The scaling relations in Eq. ([5]) define the characteristic 
time t tr when the transitions become suppressed. Indeed 
the value of the integral over time in Eq. is defined by 
the time-scale at which the phase factor starts to oscillate 
rapidly. Thus, we obtain an estimate 



a/v^Sa 2 ^) 1 /^. 



(6) 



Note that for K < 2 the loading time tf = 1/a is much 
longer than the transition time t tr provided that ha <C 
m s (Vf), which justifies the extension of the upper limit 
of integration over time in Eq.([3]) to infinity. As one 
approaches the K-T point K = 2, the loading time tf 
that obeys the relation becomes exponentially large. At 
K > 2, the system is gapless m s = 0, and the scaling in 
Eq. ([3]) is non-universal. Putting these results together we 
obtain the scaling (Q} of the excitation probability with 
the ramping rate. 
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FIG. 3: Particle-hole excitations in the Tonks-Girardeau limit 
K = 1 in real and momentum spaces in the process of loading 
of the optical lattice. The band-gap opens linearly with time. 



Similarly we can analyze the problem of merging two 
condensates with tunneling amplitude changing linearly 
in time, J(t) = Jfat, with S = \JfOi\ being the ramping 
rate. The system is described by the Hamiltonian @ 
where is the relative phase of the condensates, with 
the substitutions: V — > J and K — > 1/(2K). This map- 
ping brings the system into the attractive regime of the 
SG-model, where the low-energy excitations consist of 
breathers. The energy gap in the spectrum of breathers 
is proportional to the soliton mass, which immediately 
leads to the same rescaling for momentum and time as in 
Eqs. ([5]) with the corresponding substitution of K. The 
number of excitations scales with the ramping rate as 



2K/(6K-l), K > 1/4. (7) 



The exponent rj' decreases as one goes from strongly- 
interacting to weakly-interacting case, indicating that the 
effects of non-adiabaticity become stronger in this limit. 

In a similar fashion one can analyze other observables, 
for example, energy added to the system in a cyclic pro- 
cess in which the lattice is adiabatically loaded and un- 
loaded (merging and splitting of two ID condensates). 
Then the spectrum in the final state is gapless and the 
heating can be found as Q = J2 P e p n pi where e p is the 
excitation energy and n p is the occupation of the en- 
ergy mode. Using similar scaling arguments it can be 
shown that Q behaves as a square of the number of ex- 
citations, Q ~ (5 2 /( 3_/f ) for the problem of loading, and 
Q ~ § iK /( &K - x ) for the problem of condensates merging. 

We illustrate our main result (JlJ by solving exactly 
the TG limit K = 1. In this case, the hard-core bosons 
are mapped to free spinless fermions [3(|. The effect 
of the periodic potential is reflected in the momentum 
space through the presence of a gap between the en- 
ergy bands. In the case of a unit lattice filling, the gap 
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opens linearly in time between the first and second bands 
(see Fig. EH). Our calculation of the excitation proba- 
bility 31 1 ob tained by mapping onto the Landau-Zener 
problem [32j leads to the exact number of excitations 



= G 



n 2 h 3/2 



AV2 



(8) 



where the ramping rate S = \V/a\, and G = 0.374 is 
an exact numerical factor. Within the perturbative ap- 
proach ([3]), we find the same parametric dependence, ex- 
cept for a slightly bigger numerical factor G' = 0.44. The 
constant G' is also reproduced from the form-factor ap- 
proach using Eqs. and (g]) at K = 1. 

Another exactly solvable limit appears when [3 <C 1 in 
Eq. @. Expanding the cosine term we obtain quadratic 
Hamiltonian with a time-dependent mass term. As a 
result, we find n ox ~ log (5 1 / 3 L) confirming scal- 
ings (TT|) and 0, in the two limits K < 1 and K > 1, 
respectively. An additional factor depending on the sys- 
tem size L indicates breakdown of the adiabatic approx- 
imation 

Our results suggest an interesting observation related 
to the heating in the process of loading bosons initially 
prepared in the gapped phase at K < 2. Instead of load- 
ing bosons into the lattice directly in the gapped phase, 
to minimize the heating, one can choose a different path. 
For example, since a small lattice commensurate with the 
boson density is irrelevant perturbation at K > 2, in the 
process of loading it can be preferable first to decrease 
the interaction to values K < 2, then load the bosons 
into the lattice in the gapless phase, and afterwards adi- 
abatically change the coupling to its initial value bringing 
the system into the gapped phase. 

We estimate the number of excitations on loading 
bosons into an optical lattice by using the exact result (|8|) 
n ox = Q.Q^ayJVf/JtfEr), where Vf is the final height of 
the optical lattice and E r = (Hit) 2 / (2ma 2 ) is the recoil 
energy (the lattice spacing a is in fim and the loading 
time tf = 1/a in ms). For typical experimental pa- 
rameters [3] a = 0.425 (im and Vf/E r = 20, we obtain 
n cx = 0.19 for tf = 1 ms, and n ex = 0.06 for tf = 10 ms. 
While this number can be relatively small in a single 
ramp process, one can alternatively probe the scaling (fT]) 
by cyclically loading and unloading bosons in the lattice. 

In conclusion, we analyzed zero-temperature adiabatic 
response of one-dimensional Bose gases in two compli- 
mentary settings of loading a gas into an optical lattice 
and merging of two gases. The number of excitations 
created in these processes scales as a power-law of the 
ramping rate with the exponent strongly depending on 
the interactions in the system (see Eqs. ([T]), 0). Our 
results can be used to probe interactions in the regime 
where conventional linear response analysis fails. 
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